Interplay between Kondo and Andreev-Josephson effects in a quantum dot coupled to 

one normal and two superconducting leads 
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We study low-energy transport through a quantum dot coupled to one normal and two supercon- 
ducting (SC) leads in a junction of Y-shape. In this geometry a crossover between Kondo dominated 
and Cooper-pairing dominated states occurs by tuning the parameters such as the quantized energy 
level td of the dot and the Josephson phase <j>, which induces a supercurrent flowing between the 
two SC leads through the dot. Furthermore, Andreev scattering takes place at the interface between 
the dot and normal lead. The low-lying energy states of this system can be described by a local 
Fermi-liquid theory for interacting Bogoliubov particles. In a description based on an Anderson 
impurity model we calculate transport coefficients, renormalized parameters and spectral function, 
using Wilson's numerical renormalization group (NRG) approach, in the limit of large SC gap. Our 
results demonstrate how the Andreev resonance level approaches the Fermi level in the crossover 
region between Cooper-pairing singlet state and strong coupling situation as or cj> are varied. The 
strong coupling situation shows a Kondo effect with a significantly renormalized resonance width. 
The crossover is smeared when the coupling between the dot and normal lead is large. Furthermore, 
asymmetry in the Josephson junction suppresses the cancellations of the SC proximity for finite <\>, 
and it favors the SC singlet state rather than the Kondo singlet. 

PACS numbers: 73.63.Kv, 74.45.+C, 72.15.Qm 



-a 

o 
o 



> 
o 

(N 
m 

d 

(N 



X 



I. INTRODUCTION 

The Kondo effect in superconducting (SC) materials 
has been one of the major topics in condensed matter 
physics over forty years. The energy gap Asc of a su- 
perconductor suppresses the conduction-electron screen- 
ing at low temperatures below the Kondo temperature 
Tk- The competition of these two effects causes a quan- 
tum phase transition (QPT) between a magnetic-doublet 
and nonmagnetic-singlet ground states, which emerge at 
Asc > Tk and A sc < Tk, respectively^- 

The QPT has also been studied intensively for quan- 
tum dots^ and experiments have been carried out for 
carbon nanotube and semiconductor quantum dots* 1 ^— 
One of the merits of quantum dots is a high tunability, 
and various types of the configurations can be examined. 
For instance, for a quantum dot (QD) embedded between 
two superconducting leads in a SC/QD/SC configura- 
tion, the competition between the Kondo and Joseph- 
son effects has also been expected to occur f 1 ^— Further- 
more, an interplay between Andreev scattering and the 
Kondo effect has been studied experimentall y 28 ! 29 and 
theoretically^2r— for a QD connected to a normal-metal 
(N) lead and a SC lead in a SC/QD/N configuration. 

An interesting extension is a junction of Y-shape, at 
which a single QD is coupled to one normal and two su- 
perconducting (SC) leads as shown in Fig.[TJ This system 
has been studied by Pala, Governale and Konig, using 
a real-time diagrammatic approach based on a pertur- 
bation expansion with respect to the tunneling matrix 
elements ! 44 ! 45 Their calculations reveal precise features 



of the Josephson current and Andreev bound states for 
finite Coulomb interaction at temperatures higher than 
Tk- The approximation they used, however, is not ap- 
plicable at low temperatures T < Tk, where the com- 
petition between superconductivity and the Kondo ef- 
fect occurs. In the Y-junction, the conduction electrons 
from the normal lead can screen the local moment in- 
duced in the QD. This changes the sharp transition be- 
tween the magnetic and non-magnetic ground states into 
a continuous crossover between a Kondo singlet and a 
Cooper-pairing singlet. Furthermore, the Andreev scat- 
tering, which takes place between the QD and normal 
lead, can be controlled through the phase difference <j> 
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FIG. 1. Anderson impurity (•) coupled to one normal and 
two superconducting leads: and U are the level position 
and Coulomb interaction. T„ = irpv 2 . with p the density 
of states of the leads, and v„ the tunneling matrix element 
(u = L, R, N). The complex SC gap A L/R = \A L/R \ e ie V« 
causes the Josephson current for finite <f> = 6 R — 6 L . 
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between the order parameters of two SC leads. This is 
because the SC proximity on the QD depends sensitively 
on the properties of the junctions, and thus on 0. Con- 
versely, the Josephson current flowing between the two 
SC leads is affected by the Andreev scattering of the con- 
duction electrons from the normal lead. 

The purpose of the present work is to study these in- 
terplays of the Kondo, Andreev, and Josephson effects 
which can be observed for the QD embedded in this three 
terminal system4£ To this end, we explore a wide re- 
gion of the parameter space of this Y-junction, varying 
the position of a quantized energy level e d of the QD 
modeled with an Anderson impurity, and also examine 
how an asymmetry of the Josephson junction affects the 
low-temperature properties. Specifically, we focus on the 
crossover between the ground states which can be clas- 
sified into a Kondo singlet and a local Cooper-pairing 
singlet according to the fixed points of Wilson's numeri- 
cal renormalization group (NRG) i 47 ' 48 

For strongly correlated systems the Coulomb interac- 
tion U is larger than the other energy scales, and for 
such cases the critical behavior near the QPT is scaled 
by a single parameter Agc/Tif- in systems with QDs, 
however, some of the parameters can be tuned experi- 
mentally, and U is not always the largest energy scale. 
Therefore, the ground-state properties depend on the 
other parameters, such as e d , U, and the hybridizations 
between the dot and leads. Specifically, for small inter- 
actions U < Age, the SC pair correlations can penetrate 
into the QD and create a local Cooper pair, consisting 
of a linear combination of an empty state and a doubly 
occupied state. The essential physics of the local Cooper 
pairing can be deduced from the fixed point of the NRG 
in the limit of Asc — > oo, where the coherence becomes 
of the order of the lattice constant. We consider in detail 
this large SC gap limit in the present work. 

The low-lying energy states of this Y-junction can be 
described by a local Fermi liquid of the interacting Bo- 
goliubov particles. This is because the normal lead, cou- 
pled to the QD, has a continuous energy spectrum at the 
Fermi level. Furthermore, the local SC pair potential A d 
is induced in the QD by the proximity effect. This A d 
also plays a central role, and it varies sensitively with e d 
and 4>. For instance, the conductance due to Andreev 
scattering can be expressed at T = in terms of phase 
shift 5 for the renormalized quasiparticlcs and the angle 
Qb of the Bogoliubov rotation determined by A^, and 
is enhanced at the crossover region between the Kondo 
singlet and local-Cooper pairing singlet4i~— We also cal- 
culate the renormalized parameters for the interacting 
Bogoliubov particles. The results of the renormalization 
factor z, the Wilson ratio R, and the renormalized An- 
dreev level Ea that corresponds to quasiparticlc peak 
position appearing in the spectral function of the QD, 
provide us with sufficient information to understand the 
ground-state properties of the system thoroughly. 

This paper is organized as follows. In Sec. [HI we in- 
troduce the single impurity Anderson model for the Y- 



junction, and provide some examples which capture typ- 
ical behavior near the QPT occurring in a SC/QD/SC 
junction with a finite SC gap. In Sec. IIII1 we give a 
local-Fermi liquid description for the interacting Bogoli- 
ubov particles in the large SC gap limit, and present the 
expressions of the correlation functions in terms of the 
renormalized parameters. Then in Sec. IIV1 we provide 
the NRG results for the spectral function, transport co- 
efficient, and renormalized parameters. A summary is 
given in Sec. |Vj 



II. MODEL AND FUNDAMENTAL ASPECTS 
A. Model 

We start with the Anderson impurity model for a sin- 
gle quantum dot coupled to one normal (N) and two 
superconducting (SC) leads, 

H = H d + H »+ J2 H t,» + H sc ■ (1) 

v=N,L,R v=N,L,R 

The explicit form of each part is given by, 
H d = & (n d - 1) + y (n d - l) 2 

= E e k C l,ka C uM> 

H sc = E( a ^ c Lt c I-h + h - c -) 

k 

+ E( A ^ c Lt c k- fe i+ H - c -)- ( 2 ) 

Here, Crf = e d + U/2, and U is the Coulomb interac- 
tion. The operator d\ creates an electron with energy 
e d and spin a at the dot, and n d = d\d a . The op- 
erator cj, ka creates an electron with the energy ek in 
the leads v (= N,L,R). The couplings between the 
dot and leads are described by the tunneling matrix el- 
ements v u , and a linear combination of the conduction 
VVa = J2k c v.ka/Nv with 7V„ the number of the states 
in each lead. We assume that the density of states 
p(e) = J2k S ( e ~ e k)/N v and r„(e) = irvl p(e) arc con- 
stants independent of the frequency e at |e| < D, where 
D is the half band- width for the leads. The complex s- 
wave BCS gap, A L/R = \& L/R \e l0L ' R for the SC leads on 
the left (L) and right (R) induces the Josephson current 
when the phase difference (j> = 6 R — 9 L is finite. In this 
three terminal system, the current J v flowing from the 
dot to the lead v is given by 
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FIG. 2. (Color online) NRG results for the ground-state 
phase diagram of the Anderson impurity connected to two 
SC leads for several value of <f> = 8 R — 6 L in the electron-hole 
symmetric case 6d = — (7/2. In the upper and lower sides 
of each boundary the ground state is a nonmagnetic singlet 
and magnetic doublet, respectively. The Josephson junction is 
assumed to be symmetric T L = T R (= Ts/2) and |A^| = \A R \ 
(= A). The normal lead is not connected = in this case. 



Here, — e denotes the electron charge with e > 0. 

The Hamiltonian H contains a number of parameter 
regimes to be explored. We mainly consider the case 
where the couplings and the amplitude of the SC gaps 
are symmetric: T L = T R (= Ts/2) and |A L | = \A R \ 
(= A), for simplicity. The asymmetry in the Josephson 
coupling T L =/= T R is also examined in the last part in 
Sec. UvTJl 
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FIG. 3. (Color online) NRG results for the Josephson current 
and (d^d±) for the Anderson impurity connected to two SC 
leads are plotted as functions of <f> for several value of U in the 
electron-hole symmetric case = — U/2. The other param- 
eters are chosen such that F L — T R (= Ts/2), \A L \ = \A R \ 
(= A), T s = 2. OA, 6 R = -6 L (= 0/2), and Y N = 0. In this 
case the critical current is given by Jc = eA/h. 

correlation (d^d^) in the impurity site as functions of 
<f> for several values of U. These ground-state averages 
vary discontinuously at the QPT, and take small nega- 
tive values in the magnetic doublet ground state for fi- 
nite SC gaps A. These two expectation values in the 
doublet state are determined by the quasiparticle excita- 
tions above the SC energy gap |w| > A, and vanish in the 
limit of A — > oo. The small negative values in the doublet 
state can be explained, for instance, using the perturba- 
tion expansion with respect to 1/A from the large gap 
limit. 



III. LARGE SC GAP LIMIT 



B. QD connected two SC leads (I 1 ^ = 0) 

Before discussing the three terminal case, we first of all 
consider a simpler case with r^r = 0, where the normal 
lead is disconnected and the QD is coupled only to the 
two SC leads, in order to review some typical features 
of the competition between the Kondo and Josephson 
effects. -^^i^r— In this case the QPT occurs as a level 
crossing of the lowest two energy states of H, and thus 
the expectation value for the Josephson current and that 
for the order correlation functions show a discontinuous 
jump at the critical point. 

In Fig. the NRG results for the phase diagram of 
the ground state in the electron-hole symmetric case 
e d = —U/2 is plotted in a U /A vs Ts/ A plane for several 
values of <f>. The upper (lower) side of each boundary cor- 
respond to the parameter region where the ground state 
is a non-magnetic singlet (magnetic doublet). These re- 
sults clearly show that the magnetic-doublet region, ap- 
pearing for large U/A or small Ts/A, expands as </> in- 
creases. Therefore, the phase difference <j> between the 
two SC order parameters tends to suppresses the Kondo 
screening; 22 ' 23 ' 25 ' 26 which in this case is carried out by 
the quasiparticle excitations above the SC energy gap. 

Figure [3] shows the Josephson current J and the SC 



We consider the large gap limit, |A^/^| — > oo, in the 
following since important features of the interplays be- 
tween the Kondo effect and superconductivity in the Y- 
junction can be observed in this case although the quasi- 
particle excitations to the continuum-energy region above 
the SC gap have been projected out. For instance, the 
Andrccv resonance state emerging inside the SC gap re- 
mains near the Fermi level, and thus the essential physics 
of the low-energy transport can be extracted from this 
case. Specifically, this limit describes reasonably the situ- 
ation where the gap is much grater than the other energy 
scales, namely \A L/R \ > max(r L , T R , T N , U, \e d \). 

In the limit of |A^/^| — > oo, the Hamiltonian H can 
be mapped exactly onto a single-channel model ; 10 ' 41 ' 49 

~H-cS = H-ds + Hd + Hn + Ht,n , (4) 
n dS = A d d\d\ + A* d d l d t , (5) 
A d = T R e l6 * + T L e l6 t = \A d \ . (6) 

Thus, the SC proximity effect becomes static in this case, 
and can be described by an additional term Tids with the 
pair potential A^ penetrating into the QD. This term 
emerges as the Cooper pairs can be transferred between 
the dot and the SC leads even for large SC gaps whereas 
the unpaired quasiparticles cannot. The amplitude of 
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Ad depends on the Josephson phase, and decreases as 
increases, 



Then, in terms of the Bogoliubov particles 7^ , the effec- 
tive Hamiltonian takes the form 



|A d | 



yjl- To sin 2 (0/2) 
+ r L , To = 



r L ) a 



(7) 

(8) 



Specifically for the symmetric coupling Tr = Tl, the 
transmission probability in the normal-state case takes 
the value % = 1, and the amplitude is given simply by 

|A d | = T 5 COS (0/2) for -7T < < 7T. 



^ C ff = (n 7 -1 - 1) + — (n 7 _i - 

OO 

+ E E*;(j]+i</V + ILc < 

j=-l a 



If 



(15) 



Here, n. 
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E, 



. Tyo-Tjo-, and e<f = Ea-U/2 corresponds 



to a bare impurity level for the Bogoliubov particles. This 
representation of the Hamiltonian clearly shows that the 
total number of Bogoliubov particles, 



A. Bogoliubov particles 



The effective Hamiltonian H cff can be transformed into 
an asymmetric Anderson model for the Bogoliubov par- 
ticles, the total number of which is conservedi 25 i 41 ' 42 In 
order to carry this out, we rewrite H. ff such that 



H 



off 



d \> d l 



U Ad 
A3 -U 



\ (n d I) 2 



E E t l [fj+lafjo + fj<rfj+U 



(9) 



j = -l a 



Here, the summation over j describes the Ht,n+Hn part 
with = d a , f 0o . = ip Nt<T , and i_i = v N . The opera- 
tors fj a for j > correspond to a Wannier basis set for 
the conduction band. The explicit expression for /-^ and 
^ can be generated successively from the initial operator 
f 0cr and the energy spectrum of the conduction band, 
carrying out the Householder transformation^ There- 
fore, no approximation has been made to obtain Eq. Q 
from Eq. (jf). 

The effective Hamiltonian H eff has a global U(l) sym- 
metry in the Nambu pseudo-spin space along the direc- 
tion n oc (| Ad I cos (2d, — A^ sinOd, £d)- Thus, one can 
choose this direction n to be the quantization axis, car- 
rying out a pseudo-spinor rotation, 
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E A = je d +\^\ 2 - 



(10) 

(11) 

(12) 

(13) 
(14) 
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(16) 



is conserved in the large gap limit as a result of the global 
U(l) symmetry. Furthermore, the Friedel sum rule holds 



(n 7 _i) = -5 

7T 



(17) 
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where 6 is the phase shift of the Bogoliubov particles 
Note that even in the case of £d = where T-L eS in the 
form of Eq. (j9]) has an electron-hole symmetry, the Bo- 
goliubov particles do not have the particle-hole symmetry 
in the sense that Ea 7^ 0, and thus (n 7) _i) ^ 1, as long 
as I A d I is finite. 

At low energies the interacting Bogoliubov particles, 
described by Eq. (TT5|). show Fermi-liquid behavior that 
is characterized by the renormalized parameters: 



8 = cot 




= 1 - 



as(w) 



N 



zT 



E A = z[E A 



duj 

-E(0) 



(7 =z 2 r n;it (0,0;0,0). 



(18) 

(19) 
(20) 



Specifically, the Kondo energy scale can be deduced from 
the renormalization factor as Tk = ttTn/^- Further- 
more, from the residual interaction U between the quasi- 
particles, the Wilson ratio R can be deduced through 



R 



U 



sin 5 



(21) 



N 



We calculate these renormalized parameters with the 
NRG through the convergence of the finite-size energy 
spectrum near the fixed point41£2 In Eq. ([Tg ]) -p0 ]) . E(w) 
and r^i.^(wi, 0J2] are the self-energy and vertex 

function, respectively, for the Bogoliubov particles, the 
retarded Green's function for which is defined by 



1 I dte 1 ' 
'0 



({7- V (i),7V}} ■ ( 22 ) 



Here, the spin suffix a is suppressed on the left-hand 
side because the Green's function for a =~\ and that for 
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4- are identical due to the SU(2) symmetry for the real 
spin. The retarded Green's function for the electrons 
on the dot can be deduced from G^(lu) via the inverse 
Bogoliubov transform, 



G M (w)= -i / dte^({d a (t),dt}) 
Jo 



sm 



6b 
2 ' 



(23) 



B. Conductance and Current 

The low-energy transport, deduced from Eq. (fTS"]) . can 
also be described by the local Fermi liquid theory. At 
T = 0, the occupation number of the electrons (rid) and 
the SC pair correlation (d^d^) in the QD are determined 
by the occupation of the Bogoliubov particles (n 7 — i) = 

H 



26/ ir defined in Eq. (JT5J), and the Bogoliubov angle 9s, 
(n d ) - 1 = (<?i 7 ,-i) - l) cos Q B , (24) 
<«W = \ (K_i> - l) e ie * sin e B . (25) 

Note that cos9 B = £d/E A , e l6d smO B = A d /EA, and 
the phase of (d^d^) is given by 0d, which coincides with 
the phase of the local gap A^. The occupation number 
of the electron (rid) — 1 and (d^d^) correspond to the 
z component and the projection on the x — y plane of 
the local pseudo-spin moment induced on the impurity 
siteiSi Specifically, in the electron-hole symmetric case 
£d = 0, the local level is given by Ea = |Ad| and thus 
the Bogoliubov angle is locked at 9 B = tt/2. 

The dc conductance g NS due to the Andreev scatter- 
ing between the dot and the normal lead can also be 
expressed, at T = 0, in terms of the phase shift 5 and 
Bogoliubov angle Qsr^ 



9ns = 4r 



N 



{G r dd (uj = 0)} 12 



4e 2 

—— sin 2 9 B sin 2 26. 



(26) 



Here, {G r dd (uj)}i2 = ((d^,;d^)) u is the off-diagonal 
(anomalous) component of the retarded Green's function 
in the Nambu pseudo-spin formalism, the corresponding 
Matsubara function of which is defined in Appendix [Al 
Therefore, the zero-temperature conductance g NS is de- 
termined by the value at the Fermi level u> = 0. The 
argument 26 in Eq. (|26j) appears as a difference between 
the phase shift for the quasi-particle +6 and that for the 
quasi-holc —6. 

In the large gap limit, the Josephson current flowing 
through the dot can also be expressed in terms of the 
Bogoliubov angle 9 B and the phase shift 6, or (n 7j _i), 
which also depend on (f> [see Appendix |A"] . 
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FIG. 4. (Color online) Ground-state phase diagram in the 
limit of Fn —> and |A L / fl | — > oo is plotted in e<i vs F s plane 
for several values of <f>. The couplings between the QD and 
SC leads is chosen to be symmetric F L — F R (= r g /2). In 
this case the local SC gap is given by = r s cos<^/2 for 
\4>\ < 7T. The ground state is a singlet due to the local Cooper 
pairing outside the semi ellipse of sj(e d + C//2) 2 + |A d | 2 = 
(7/2. Inside this semi ellipse the ground state is a magnetic 
doublet, which for finite Ym changes to a Kondo singlet due 
to the conduction-electron screening of the local moment. 



Here, Jc = eTs/h is the critical current. 



C. Kondo singlet vs Local Cooper pairing 

The ground state of the asymmetric Anderson model, 
given in Eq . (fl5|) , can be classified according to the fixed 
points of the NRGj 47 ' 48 Among them the strong- coupling 
fixed point describes the Kondo singlet, for which the im- 
purity site is singly occupied. The frozen-impurity fixed 
point describes a different situation, where the impurity 
level is far away from the Fermi level and the impurity 
site becomes empty or doubly occupied. In our case, the 
frozen- impurity fixed point is defined with respect to Eq. 
([To]) for the Bogoliubov particles, and thus this describes 
a singlet state by a local Cooper pairing that consists of 
a linear combination of the empty and doubly occupied 
impurity states of the electrons represented in Eq. (|9|). 
We refer to this fixed point as local Cooper pairing in the 
following. 

The ground state of ~H c g varies continuously between 
the Kondo singlet and the local-Cooper pairing, depend- 
ing on the Hamiltonian parameters Ea, U, and T^- A 
rough sketch of the ground-state phase diagram can be 
obtained quickly from that in the atomic limit T N — > 0, 
where the normal lead is disconnected. In this limit the 
dot is occupied by a single Bogoliubov particle with spin 
1/2 for Ea < U/2 whereas the dot is empty with no Bo- 
goliubov particle for Ea > U/2. Thus, the phase bound- 
ary is given by a simple semi ellipse of + \Ad\ 2 = 
U/2, which is plotted in Fig. |4] for several values of <j>. 
The ground state is a singlet due to the local Cooper 
pairing outside the semi ellipse whereas inside the semi- 
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elliptic boundary the ground state is a doublet and thus 
the local moment arises in this limit of Tjv — > 0. How- 
ever, the local moment is screened when the normal lead 
is connected F^r 7^ 0, and then the ground state inside 
the semi ellipse changes to the Kondo singlet. The cou- 
pling to the normal lead also changes the sharp transition 
at the boundary to a continuous crossover between the 
local-Cooper-pairing singlet and the Kondo singlet4i 

There are further quantitative corrections when the 
SC gap A is finite. This was also examined for the IV = 
case, using the NRG&i The results showed that the 
region of the magnetic doublet state becomes small as 
A decreases. This is because also the excited states in 
the SC leads above the gap \u\ > A contribute to the 
screening of the local moment for finite A. 



IV. NRG RESULTS 

In this section, we provide the NRG results for the 
ground-state properties of the Y-junction in the large gap 
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A. Spectral function 

We first of all discuss the impurity spectral function for 
electrons Add = (~ l/ 71 ") FnG^ that can be deduced from 
the one for the Bogoliubov particles A 1 = (— I /w) ImG 7 , 

A dd {u) = A» cos 2 ^ + A,(-u) sin 2 % . (28) 



Specifically, the spectrum is symmetric A dd (— uS) = 
A dd (uj) in the electron- hole symmetric case, where &b = 
7r/2. Note that the single-electron spectrum A dd (uj) con- 
sists of a superposition of a single-Bogoliubov-particle 
part A^(uj) and a singlc-Bogoliubov-hole part A^{— lu). 
This can be deduced from the Lehmann representa- 
tion that is expressed in terms of the matrix element 
(to, iVj|dl|GS, N b ) between the ground state GS, N b ) and 
an excited state |to, of ~H e e, where N b is an eigenvalue 
of the total number of the Bogoliubov particles 7V 7 de- 
fined in Eq. (|16p . This matrix element can be finite not 
only for the single-particle excitations with N b = N b + 1 
but also single-hole excitations with = N b — 1 of the 
Bogoliubov particles because the electron d j can be de- 
composed into a superposition of the annihilation 7^ ^ 
and creation 7_ 1 , of the Bogoliubov particles. 

The low-energy spectral weight is dominated by the 
rcnormalizcd Andrcev resonances that appear in A dd {u) 
as a pair of the quasiparticle peaks, 



T N cos 2 % 



r 7v sin 2 % 



oj - E A 



E A 




(29) 



FIG. 5. (Color online) NRG results for the Spectral func- 
tion A dd {uj) = (-l/ir)lmG ddi!r (u), in the large gap limit 
l^i/fll — > 00 a t <fi — 0, is plotted for several value of U choos- 
ing the couplings to the SC leads such that (upper panel) 
rs/Tjv = 5.0, and (lower panel) r s /riv = 1.0. The other 
parameters are taken to be Fl = Tr (= Ts/2) assuming the 
electron-hole symmetric e d = — U/2. In the present case the 
local SC gap in the impurity site is given by Ad = Tg, and 

E A = r s . 



The peak position iE A , width Tn, and renormalization 
factor z vary as the Coulomb repulsion increases from 
U = 0, for which we have E A = E A , Tjv = IV, and 

2=1. 

The high-energy part of the spectral weight away from 
the Fermi level can be inferred from the one in the atomic 
limit r^y = 0. For weak repulsions E A > U/2, the 
ground state is a singlet, and A dd (uj) has two peaks at 
ui = ±(E A — U/2). This is because in this case A 7 (u) has 
a single peak at u = E A — U/2, which moves towards the 
Fermi level as U increases. On the other hand, for strong 
repulsions E A < U/2 in the T N —> limit, the ground 
state is a magnetic doublet, and then A dd (to) has four 
peaks emerging at u = ±£xjp anc ^ = ' =£ lw- These peaks are 
caused by the excitations to the upper e UP = U/2 + E A 
and lower e LW = —{U/2 — E A ) atomic peaks defined 
with respect to the Bogoliubov particles, and each of 
the two final states consists of a linear combination of 
an empty and doubly occupied states of the original elec- 
trons. These peaks correspond to the original Andrcev 
typically found in situations of a QD coupled only to a 
superconductor, and they are delta-functions in the limit 
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= 0. For finite superconducting gap their position 
and the occurrence of the ground state transition change 
quantitatively.— In this situation in the doublet phase it 
is possible that the higher excitations ±£up are n0 ^ f° un d 
within the gap anymore. 

Figure [5] shows the spectral function for <fr = in the 
electron- hole symmetric case, where Ea = T s . In the up- 
per panel A dd (oj) in the case of relatively small T N with 
Ea = 5T N is plotted for several values of the Coulomb 
repulsion U/T N = 0, 5, 10, 15, 20. In these examples, at 
U = lOr^ the bare parameter takes the value Ea = U /2, 
and thus the dash-dot line represents the results obtained 
at the crossover region between the Kondo and local SC 
singlet states. The pair of the renormalizcd Andreev res- 
onances at oj ~ ^E A , which correspond to the Kondo 



peak for the Bogoliubov particles described in Eq 
shift closer to the Fermi level oj = as the Coulomb re- 
pulsion U increases from to 2Ea- Then, for U > IE a 
five peaks emerge as seen the curves for U = 15T N and 
20T N . Among them, the central peak near the Fermi 
level w = corresponds to the Kondo resonance for the 
Bogoliubov particles, which appears for A 1 (oj) at small 

positive frequency u ~ E A . However, for the electron 
spectral function A dd (oj), this peak at oj ~ E A and the 
counterpart for holes at oj ~ — E A overlap to form a sin- 
gle peak at the Fermi energy oj = 0. The other four 
peaks represent the electron and hole components of the 
excitations to the the upper and lower atomic levels of 
the Bogoliubov particles. For instance, in the curve for 

-5Fjy and the one at 



U = 20T N , the broad peak at oj 
lu ~ +15r A , correspond to upper and lower atomic peaks 
appearing in A^(ui), respectively. 

The lower panel of Fig. [5] shows the spectral function, 
obtained at Ea = F N where the bare Andreev level and 
the hybridization energy due to the coupling to the nor- 
mal lead coincide. The Coulomb interaction is chosen to 
be U = 5T N and lOr^, so that Ea < U/2 for these two 
cases. As the coupling to the normal lead T N is relatively 
large in these cases the sub peaks of the Coulomb oscil- 
lation are smeared especially in the curve for U = 5T N . 
Nevertheless, in the other curve for U = 101^, the sharp 
central peak and two sub peaks of the atomic nature are 
distinguishable. The central peak emerges as a results of 
the superposition of the two renormalized Andreev level 
at oj ~ ±Ea close to the Fermi level whereas each of the 
sub peaks at oj ~ ±4.0Fat emerge as a sum of the sub 
peaks in A 7 (oj) and that in the counterpart for the holes 

Ay(-W). 

The results which we have presented in Fig. [5] have 
been the ones obtained at = 0, where there is no phase 
difference between the two SC order parameters. As we 
have chosen the parameters such that £d = and Tl = 
Tr, the dependence of Add{oj) on <f) enters only through 
the bare Andreev level Ea = Tscos(f>/2 in this case. 
Therefore, the results capture essential features common 
to the case for finite <f>. In the electron-hole asymmetry 
case £d 7^ 0, however, Add{w) is no longer a symmetric 




FIG. 6. (Color online) Phase shift and some related ground- 
state averages are plotted vs Josephson phase cj> for several 
values of Tjv in the electron-hole symmetric case e c i = —U/2: 
(upper panel) phase shift 5 and pair correlation 2 I (d^d-j-) I , 
(lower panel) conductance g NS and Josephson current J in 
units of Jc = er s /ft. The other parameters are chosen such 
that r_L = F_r (= Ts/2), and U = 1.5Fs in the large gap limit 
\Al/r\ -> oo. 



function of u. The asymmetry in the oj dependence enters 
through Qb as it deviates from 7r/2 for £<j ^ 0. 



B. Transport properties & Fermi-liquid behavior 

In this subsection we present the NRG results for the 
ground-state properties of the Y-junction for symmetric 
coupling = (= r s /2), where To = 1- The Joseph- 
son phase is chosen such that Br = —6l (= <fi/2), which 
makes the local SC gap and pair correlation (d^d^) 
real 0d = as shown in Appendix. We consider the 
electron-hole symmetric case £d = in Sec. IIVB li and 
then discuss the gate voltage dependence varying ^ in 
Sec. IIVB 21 



1. Electron-hole symmetric case 

Figure [6] shows the results of the phase shift S, the SC 
pair correlation 2|(d,(i^}|, the conductance g N g, and the 
Josephson current J as functions of cj> for several values 
of Tat/Fs in the electron-hole symmetric case, where the 
Bogoliubov angle is locked at 8s = 7r/2 and (rid) = 1, as 
mentioned. Thus, the ground state properties are deter- 
mined by the phase shift 5, which depends on the Joseph- 
son phase through Ea = T s cos(0/2). The Coulomb 

1.5rg, and for this in- 
- 0.46tt in the T N -> 



repulsion is chosen to be U 
teraction the QPT occurs at 
limit. 

The coupling to the normal lead makes the excitation 
spectrum of the dot gapless and changes the sharp QPT 
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FIG. 7. (Color online) Renormalized parameters are plotted 
vs Josephson phase 4> for several values of Fjv in the electron- 
hole symmetric case e<j = —U /2: (upper panel) renormalized 
Andreev level ±Ea, (lower panel) renormalization factor z 
and Wilson ratio R. The other parameters are chosen such 
that r_L = r_R (= Ts /2), and U = 1.5Fs in the large gap limit 
\A L/R \ ao. 



into a continuous crossover between the two different sin- 
glet states, namely the Kondo and local SC singlets. For 
<f> < 0.467T, the ground state is the local Cooper pair- 
ing consisting of a linear combination of the empty and 
doubly occupied impurity states with small S. For large 
phase difference </> > 0.46-7T, however, the ground state is 
a Kondo singlet state, for which 6 ~ ir/2. 

The phase shift shown in the upper left panel of Fig. [5] 
can be regarded as the number of Bogoliubov particles on 
the impurity site, (n 7l _i) = 25/tt, due to the Friedel sum 
rule given in Eq. (|I7I) . Furthermore in the electron- hole 
symmetric case, the SC pair correlation defined in Eq. 
([21)]) is given simply by 2(tLdi ) = 1 — (n 7i _i). Therefore, 
the pair correlation is suppressed in the Kondo regime 
for 4> si 0.467T, as seen in the upper right panel. 

The Andreev conductance g N g, shown in the lower left 
panel of Fig. [51 also varies as a function of the Josephson 
phase <$>. At the crossover region between the Kondo 
singlet and local-Coopcr-pairing states, the conductance 
due to the Andreev scattering g NS has a sharp peak for 
small T N . The crossover behavior, however, is smeared 
as r^r increases. The conductance takes the unitary limit 
value 4e 2 /h at 5 = 7r/4 where the renormalized resonance 
widtli^and the renormalized Andreev level coincide such 
that r„ = E A . This happens in Fig. [5] at <f> ~ 0.0 for 



N 



o.sr 



s- 



The Josephson current, in the lower right panel, also 
shows the crossover behavior near (f) ~ 0.467T, and de- 
creases at <f) > 0.467T when T^r is small. The value of the 
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FIG. 8. (Color online) Phase shift and some related ground- 
state averages are plotted vs ed/U for <j) = 0.6-7T for several 
values of Fjv: (upper panel) phase shift 8 and electron oc- 
cupation number (riti), (lower panel) conductance g NS and 
Josephson current J in units of Jc = eT s / h. The parameters 
are chosen to be T L = T R (= r s /2), and U = 1.5r s . 



current approaches zero in the Kondo-singlet region since 
in this case the large gap limit has been taken. For finite 
SC gaps, however, a weak current will flow in the oppo- 
site direction as seen in Fig. [3] for the magnetic-doublet 
state. Similarly, for finite SC gaps, the SC pair corre- 
lation (d^dj) will also change the sign at crossover, and 
has a small negative value in the magnetic ground state. 

We have also deduced the renormalized parameters for 
the local Fermi liquid of the interacting Bogoliubov par- 
ticles from the convergence of the finite-size energy spec- 
trum in the successive NRG steps . 47 ' 52 Figure [7] shows 
the results for the renormalized Andreev level ±E A , the 
wavefunction renormalization factor z, and the Wilson 
ratio R. We see for small r^v (= 0.05Fs) that the param- 
eters are significantly renormalized in the Kondo singlet 
state for <fi > 0.467T, where z <Si 1.0, R ~ 2.0, and the pair 
of the renormalized Andreev peaks ±£U lie close to the 
Fermi level io = 0. This indicates that the Bogoliubov 
particles are strongly correlated in the Kondo regime. In 
contrast, in the local Cooper-pairing state for <f> < 0.467T, 
the parameters are not renormalized so much z ~ 1.0, 
R ~ 1.1, and the renormalized Andreev peaks ±E A sit- 
uate away from the Fermi level. As the coupling r^r 
between the QD and normal lead becomes large, these 
two singlet states become indistinguishable, as those for 
T N = 2.5T S seen in Fig. □ 



2. Electron-hole asymmetric case 

We next consider the electron-hole asymmetric case, in 
which the Bogoliubov angle &b = cot - 1 Ad|) devi- 
ates from 7r/2 as £<j varies. Furthermore, the phase shift 
6 also varies as a a function of £d since the bare position 



FIG. 9. (Color online) Phase shift and some related ground- 
state averages are plotted vs e<j/t/ for cj> = 0.3-7T for several 
values of Fat: (upper panel) phase shift 5 and electron oc- 
cupation number (rid), (lower panel) conductance g NS and 
Josephson current J in units of Jc = eT s /H. The parameters 
are chosen to be Fl = Fjj (= Ts/2), and U = 1.5Fs- 



FIG. 10. (Color online) Renormalized parameters are plotted 
vs ed/U for several values of <j> for relatively small coupling 
Fn — 0.05Fs between the dot and normal lead: (upper panel) 
renormalized Andreev level ±Ea and phase shift S, (lower 
panel) renormalization factor z and Wilson ratio R. The other 
parameters are chosen such that Tl = Tr (= Ts/2), and 
U = 1.5F S . 



Ea — U/2 of the Andreev level depends on Through 
these changes of the phase parameters 0b and 5, the 
ground-state properties of this Y-junction depend on the 
gate voltage e^. 

In this subsection we examine the td dependence for 
several values of the Josephson phase <f) = 0.37T, 0.467T, 
and 0.67T, choosing the Coulomb repulsion to be [/ = 
1.5Ts as that in the half- filled case discussed in the above. 
The phase boundary between the singlet and doublet 
ground states moves in the T s vs ed plane as 4> increases 
as shown in Fig. |4] In our parameter set Ts/U = 0.666, 
and the QPT occurs for <fi = 0.6-7T when a varies in the 
range —1.0 < td/U < 0.0, whereas the level crossing docs 
not occur for cf> = 0.3tt. A marginal situation is realized 
for <f> ~ 0.467T, in this case the system approaches very 
closely to the phase boundary near the symmetric point 
td * -U/2. 

In Fig. [51 the correlation functions for <f> = 0.6n are 
plotted vs ed/U. The ground state changes discontinu- 
ous^ in the limit of T N = 0, at e d ^ -0.82[/and -0.18/7. 
The sharp transition becomes a continuous crossover for 
finite T N , and at -0.82C7 < e d < -0.18Z7 the ground 
state is a Kondo singlet consisting of the strong corre- 
lated Bogoliubov particles, and the electron filling is al- 
most (rid) — 1.0, as shown in the upper right panel. In 
the upper left panel, we can see that in this region of 
ed the occupation number of the Bogoliubov particles at 
the impurity level, (n 7 _ x ) = 25/ix, decreases as T N in- 
creases from zero to an intermediate value T N < 0.5r s . 
Then, (n x ) increases as seen in the results obtained for 
= 2.5r s , and approaches 1.0 in the limit of large T N 
where this coupling dominates all the other effects. The 
lower right panel of Fig. [8] shows the Josephson current. 
This supercurrent, flowing between the two SC leads, is 



suppressed due to the electron correlation in the Kondo 
regime at -0.82U < e d < -0.18C/. Out side of this re- 
gion the ground state is characterized by the local Cooper 
pairing, and the current is not so suppressed although the 
coupling to the normal lead T N smears the structure due 
to the QPT. We also see in the lower left panel that the 
conductance g NS between the normal lead and the QD 
shows sharp peak at the transient region of the crossover 
for small T N . The sharp conductance peak is mainly 
caused by the phase shift 8 that changes suddenly from 
to 7r/2 at the crossover region because the conductance 
is proportional to sin 2 25. The Bogoliubov angle 8s, ap- 
pearing in the expression of g NS given in Eq. (|26p , varies 
moderately and determines the peak height. 

Figure [5] shows the ed dependence of the correlation 
functions for a smaller value of the Josephson phase 
<f) = O.Stt. In this case U is not large enough to reach 
the Kondo regime over the crossover region. The ground 
state is the local Cooper-pairing state for all values of 
ed, and thus the correlation functions vary moderately as 
a function of the gate voltage ed- The conductance and 
Josephson current have maximum at the electron-hole 
symmetric point = — 0.5U. This is mainly because 
the factor sin Qb that appears in the expression of these 
correlations given in Eqs. (|26[) and (|2"T|) takes a local max- 
imum at the Bogoliubov angle of 0^ = tt/2. 

In Fig. [TU1 the results of the renormalization factor 
z, the renormalized Andreev levels ±Ea, the Wilson ra- 
tio R, and the phase shift are compared for three dif- 
ferent values for the Josephson phase cf> = O.Stt, 0.467T, 
and 0.67T. The Coulomb interaction and the hybridiza- 
tion energy scales are chosen such that U = 1.5r,5 and 
Tat = 0.05rs- We can see in the upper left panel that the 



10 



pair of renormalized Andreev levels ±Ea for </> = 0.67T lie 
very closely to the Fermi level at -0.8217 < e d < -Q.18U. 
Furthermore, in this region, the renormalization factor z 
significantly decreases and the Wilson ratio approaches 
R — > 2 owing to the strong correlations in the Kondo 
regime. It also indicates that the Kondo temperature 
Tjc = ttYn/4 and the renormalized resonance width r^r, 
defined in Eq. (|18p . become very small. Simultaneously, 
the local Bogoliubov particles in the dot approaches the 
single occupancy (n 7i _i) = 25/ir ~ 0.9 although it is less 
than 1.0 because T N is not very small in this case. Note 
that the sharp single Kondo peak that we have seen in 
Fig. E] for U > 15Tn consists of such a pair of the renor- 
malized Andreev levels, appearing in the close vicinity of 
the Fermi level. 

The ground-state properties show a marginal behav- 
ior at <j> — 0.46-7T, as shown in Fig. [TO] For instance, z 
and R take the intermediate values, and the pair of ±Ea 
becomes distinguishable as we can see in the upper left 
panel of Fig.[]I)J Then, for a smaller value of the Joseph- 
son phase cj) = 0.3-7T, the ground state is a singlet caused 
by the local Cooper pairing, as mentioned. Therefore, 
in this case the electron correlations are suppressed as 
z ~ 0.95 and R ~ 1.16 even at the electron-hole symmet- 
ric point ed = — 0.5f7. 



C. Spatial asymmetry in the junction (Yl 7^ Yr) 

So far, we have assumed that the Josephson junction is 
symmetric Yl = Yr. However, the SC proximity effects 
in real systems depend on the asymmetry in the couplings 
Tl 7^ Fr. One difference of the asymmetric junction 
from the symmetric one is that the maximum possible 
of the transmission probability is no longer one, namely 
7o < 1. Specifically, the local SC gap induced in the dot 
I Arf I , defined in Eq. , becomes larger in the asymmetric 
junction than that in the symmetric junction for which 
7o = 1. Therefore, the asymmetry in the couplings tends 
to enhance the SC proximity effects as it suppresses the 
reduction of |A<j| due to the Josephson phase <j>. In the 
following, we discuss the ground state properties of the 
asymmetric junction with Yl =/= Fr, keeping the sum 
Ys = Fr+Tl unchanged at Yn = 0.05Fs. For simplicity, 
we examine the electron-hole symmetric case ed = —U /2, 
where the bare Andreev level is given by Ea = |A<j| and 
the Bogoliubov angle is fixed at 6_b = tt/2. 

The NRG results for the asymmetric junction are plot- 
ted vs the Josephson phase </> in Fig. [IT] for U = 1.5rg. 
In each panel, the curve for Yl/Yr = 1.0 corresponds 
to the results obtained for the symmetric coupling, pre- 
sented also in Fig. [6] As mentioned in the above, the 
amplitude of the static SC gap |A^| for the asymmetric 
junction becomes larger that for the symmetric junction. 
Thus the position of Ea moves away from the Fermi en- 
ergy due to the spatial asymmetry, and this causes rather 
moderate <j> dependence of the phase shift 8 for Yl 7^ Tr, 
seen in the upper left panel of Fig. 111! Correspondingly, 
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FIG. 11. (Color online) Phase shift and some related 
ground-state averages for the asymmetric couplings Tr 7^ Yl 
are plotted vs Josephson phase <f> for several Yl/Yr keeping 
Ys = Yr + Yl unchanged: (upper panel) phase shift 5 and 
pair correlation 2 | (d±d^) | , (lower panel) conductance g NS and 
(lower right) Josephson current in units of Jc = eY s /h. The 
other parameters are chosen such that ed = —U/2, U = 1.5Fs 
and Y N = 0.05Y S . 



the SC pair correlation 2(d^d^}, shown in the upper right 
panel, increases as the asymmetry T^/Yr increases. 

The conductance due to the Andreev scattering is pro- 
portional to sin 2 25, and thus has a peak when the phase 
shift takes the value of S = 7r/4 at crossover region be- 
tween the local-Cooper pairing and Kondo singlet states. 
We see in the lower left panel of Fig. [TT] that the posi- 
tion of the conductance peak shifts towards the larger <p 
side as 5 becomes smaller with increasing Yl/Yr. This 
also shows that the asymmetric coupling favors the SC 
proximity into the dot, and enlarges the parameter region 
for the local-Cooper-pairing ground state. Furthermore, 
the crossover behavior from the local-Cooper pairing to 
the Kondo regime, seen for the Josephson current in the 
lower right panel for Yl/Yr = 1, is smeared as Yl/Yr 
increases, and the current shows a simple sinusoidal <fi 
dependence for large asymmetries. 



V. SUMMARY 

We have studied the crossover between a Kondo singlet 
and a local-Copper-pairing singlet, occurring in a quan- 
tum dot coupled to one normal and two SC leads. The 
low-energy states of the system can be described by ex- 
citations from a local Fermi-liquid ground state of inter- 
acting Bogoliubov particles. Specifically in this three ter- 
minal configuration, the renormalized parameters for the 
quasiparticles vary as functions of the Josephson phase c/>, 
and the crossover occurs at finite <p c as the phase differ- 
ence varies in the range < \<j>\ <ir. We have calculated 
the phase shift 8, the renormalization factor z, the renor- 
malized Andreev level ±Ea, and the Wilson raton R in 
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the large SC gap limit Age — > oo, using the NRG, and 
have deduced the transport properties at T = 0. 

The Bogoliubov particles are strongly renormalizcd in 
the Kondo regime while the renormalization is a minor 
effect in the local-Cooper-pairing regime which corre- 
sponds to the frozen- impurity fixed point of the NRG. 
Near the crossover between the two regimes a pair of the 
renormalizcd Andreev levels ±Ea approach the Fermi 
level, and the conductance between the dot and normal 
lead has a peak. The Joscphson current between the two 
SC leads is suppressed significantly in the Kondo regime. 

We have also presented the spectral function calculated 
with the NRG. The results demonstrate precise features 
of the original Andreev levels, which for the local-Coopcr- 
pairing state arc broadened by Tm and then rcnormal- 
ized as U increases. In the Kondo regime, the pair of 
the renormalizcd Andreev levels overlap to form a single 
peak near the Fermi level while at high energies four ad- 
ditional peaks are visible. These peaks correspond to the 
excitations to the upper and lower atomic peaks, which 
are defined with respect to the Bogoliubov particles and 
consist of a linear combination of the empty and doubly 
occupied electron states. Thus, the broadened bare An- 
dreev peaks and a low energy feature corresponding to a 
Kondo resonance can appear within the superconducting 
gap. For suitable parameters this should be observable 
experimentally in the discussed Y-shapc geometry. 
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Here, ioj is the Matsubara frequency, is the self 

energy due to the Coulomb interaction, g v (iu) is the local 
Green's function at the junction of the lead v 



9v(. ibJ ) = -Kp v 



iuj\ — A„ 



A„ 
A* 



(A3) 



1 and T3 are the unit and Pauli matrices, respectively. 

The Josephson current from the dot to the SC lead 
for v = L,R can be expressed in terms of the Green's 
function 



h p 

X 



^2 Tlc [{9A i ^ J n)T 3 - T 3 g u {iu] n )} G dd (iu) n ) 

(A4) 



In the limit of lA^I — > 00, the lead Green's function 
becomes a constant g u — > irp u A^/l A„| as the retarda- 
tion effects caused by the u> dependence in Eq. (|A3[) are 
suppressed. Then, Eq. (|A4[) can be rewritten in the form, 
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Appendix A: Current in the large gap limit 
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(A5) 
(A6) 



Note that (c^df) = \(did-f)\ e t6d , in the limit of large gap 
as the phase of (d|,dt) coincides with that of the local SC 
gap Ad = |Ad|e l6,<i as shown in Eq. ([23)1 . The current 
conservation (Jr) + {Jl) = can be confirmed explicitly 
through the identity 



We provide the expression of the Josephson current in 
the large gap limit in this appendix. To this end, we use 
the imaginary time Green's function, defined by 



G m (t) 



(T T dJr)dh (TrdJ^dJ 
(T T d\(T)d\) (T T 4(r)^) 



(Al) 



The Fourier transform of this function can be expressed 
in the form, 

^G d d(iu)} = iojI - £ d T 3 + iT N sgnul 

- Y, v 2 u9v(i") - SH. (A2) 

v=L,R 



T R sin (6 R -e d )- T L sin (6 d -0 L ) - , ( A7) 



which follows from the definition of 9d given in Eq. ([5]). 
Using Eq. (|A7[) , the current can be expressed in the form 



(A8) 



This can be rewritten further, in terms of the phase shift 
S and the Bogoliubov angle 6s , as shown in Eq. (|27p . 
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